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Method of and apparatus for error correction of multiplex signals 

Cross Reference To Related Applications 

This application is a continuation of International Application 
PCT/EP03/11193 having an international filing date of Oct. 9, 2003, and 
from which priority is claimed under all applicable sections of Title 35 of 
5 the United States Code including, but not limited to. Sections 120, 363 

and 365(c), and which in turn claims priority under 35 USC 5119 to 
German Patent Application No. 102 48 507.0-42 filed on Oct. 11, 2002. 

Technical Field 

Reducing errors in telecommunications, and more particularly in 
10 OFDM type telecommunication signals. 

Background Art 

The multi-carrier method OFDM (orthogonal frequency division 
multiplex) is becoming more widespread as time goes by. It is used for 
digital radio (digital audio broadcasting or DAB), digital television (digital 

15 video broadcasting or DVB) and for local radio networks (high 

performance local area network or Hiperlan). 

OFDM differs from single-carrier methods in that the information 
flow is transmitted not by way of a single carrier but by way of a number 
of subcarriers. A frame to be transmitted is subdivided into a number of 

20 symbols. Each symbol contains a number of data which are transmitted 
from the transmitter to the receiver distributed to various subcarriers. A 
given frequency spacing is maintained between the subcarriers. The 
orthogonality of the subcarriers is achieved by orthogonally encoded 
digital modulation. The orthogonality of the subcarriers ensures 

25 distinguishability thereof even in the case of spectral overlap. 

Not all subcarriers in an OFDM system contain user information. 
There are so-called pilot subcarriers, also referred to hereinafter as pilot 
channels, for the transmission of pilot symbols. The pilot symbols contain 
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information which is known from the outset to the receiver. In the 
transmission of a frame a plurality of pilot signals are transmitted, which 
are integrated into the frame at usually regularly spaced positions. 

Pilot subcarriers are basically used for channel estimation. In that 
5 case channel attenuation of the various subcarriers which are also referred 
to as subchannels is ascertained by means of interpolation of the received 
values of the pilot subcarriers. The standards "IEEE 802.11a" and 
u Hiperlan/2" for wireless communication in networks in the 5 GHz 
frequency range do not permit use of the pilot channels for channel 

10 estimation. The pilot channels are here provided in mutually different 
frequency ranges which involve a spacing of 4.375 MHz. Therefore, only 
channels with a maximum delay spread of 2/4.375 MHz = 0.46 \as can be 
correctly estimated by interpolation of the pilot channel information. In 
general terms however the pilot subcarriers cannot be used in the 

15 specified standards for channel estimation. 

Correct decoding of OFDM signals on the receiver end is delicately 
dependent on the efficiency of the synchronization unit of the receiver. 
The synchronization unit is responsible for the detection of incoming 
frames and for estimating and correcting possible frequency offsets. 

20 Frequency offsets produce errors which lead to incorrect decoding of the 

received signal. 

WO 01/20863 Al describes a method of and an apparatus for 
correction at the receiver end of a phase error of the received signal in the 
time domain. There is provided a phase estimation circuit having two 

25 phase locked loops. The correction of a phase error in the time domain 

has the disadvantage of a great delay for the signal. The use of phase 
locked loops for phase estimation additionally signifies a high level of 
circuitry complication and expenditure. 

Other methods provide for phase estimation and phase correction in 

30 the frequency domain. Provided for that purpose is an FFT unit which 
transforms the received signal by means of a fast Fourier transform or FFT 
into a frequency spectrum which can be broken down according to 
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subcarriers. In those methods a synchronization unit also decides about 
the start time from which the incoming signal is subjected to a FFT. 

The publication V. Mignone, A. Morello: "CD3-OFDM: A Novel 
Demodulation Scheme for Fixed and Mobile Receivers" (IEEE Transactions 
5 on Communications, Vol. 44, No. 9, September 1996, pages 1144 to 

1151) discloses a method and an apparatus in which the Fourier- 
transformed signal of each subcarrier is subjected to equalization prior to 
decoding. Equalization is effected by means of the signal received in the 
preceding cycle. For that purpose that signal is decoded, then recoded in 

10 a kind of feedback loop and used for channel estimation of the 
corresponding subcarrier. The frequency response function ascertained 
from the channel estimation is used for equalization of the current 
Fourier-transformed signal by the current signal being divided by the 
frequency response function ascertained from the preceding signal. 

15 It has been found however that, after the Fourier transform and 

equalization in the frequency domain and decoding of the symbols phase 
errors are contained in the equalized signal. That gives rise to errors in 
the decoding operation, which can result in a data communication 
crashing. 

20 Conventional methods of eliminating that residual phase shift 

involve correction of the phase error by means of complex multiplication 
with a suitable phase factor of the value 1 (phasor). Complicated circuits 
which are based on numerically controlled oscillators (NCO) are usual for 
correction of the phase error. Those circuits cause a severe delay in 

25 terms of signal decoding. 

Disclosure of the Invention 

The technical object of the invention is to provide a fast method of 
reducing a phase error in a received signal, which is caused by at least 
one error source, and a suitably designed apparatus for carrying out the 
30 method. 
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That object is attained by a method of reducing a phase error 
caused by a plurality of error sources in a signal which is present in a 
digital frequency representation in the form of a sequence of a plurality of 
digital partial signals which are associated with a number of subcarriers 
5 (k) of a carrier, with the steps which are effected for each partial signal: 

- equalization of the partial signal, 

- estimation of the phase error of the equalized partial signal, and 

- correction of the estimated phase error of the equalized partial 

signal. 

10 In the method according to the invention the equalization step 

includes a step of at least partially eliminating an accumulation of a phase 
error of the partial signal, caused by a sampling frequency error, over the 
sequence of the partial signals, such that the accumulation is negligible. 
The estimation step includes a step of detecting a plurality of 

15 predetermined pilot signals and a step of determining a phase correction 

factor on the basis of the detected pilot signals, wherein at least one 
multiplication operation is carried out solely by means of shift and adding 
operations. 

Brief Description of the Drawings 

20 The method according to the invention is based on an analysis of 

the responsible error sources for the phase errors remaining after the 
synchronization and Fourier transform steps in an OFDM receiver. In that 
respect reference is firstly made hereinafter to Figure 1 to describe the 
structure of a typical OFDM receiver in accordance with the state of the 

25 art. The essential sources of the remaining phase errors are then 

discussed with reference to Figures 2 and 3. 

Detailed Description of the Invention 

Figure 1 is a simplified block circuit diagram showing a per se 
known OFDM receiver in the form of a heterodyne radio frequency down- 
30 converter 10, hereinafter referred to for brevity as the down-converter. 
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The down-converter 10 has an antenna 12. Connected downstream 
of the antenna is a band pass filter 14 and a low-noise pre-amplifier 16 
(referred to as a low-noise amplifier or LNA). 

A first mixer stage which is connected at the input side to the 
5 output of the pre-amplifier has a frequency mixer 18 and a local oscillator 

20. The frequency of the oscillation produced by the local oscillator 20 is 
the difference of the carrier frequency F c of the oscillation received by the 
down-converter at its antenna and an intermediate frequency fi F (referred 
to as IF). The frequency mixer multiplies the filtered and pre-amplified 

10 signal by the frequency of the local oscillator 20. Due to tuning 

inaccuracies in the local oscillator 20, a frequency offset of the frequency 
of the local oscillator 20 of ±f E can occur. That is considered hereinafter 
with reference to Figure 2. 

A channel filter 22 is connected on the output side of the local 

15 oscillator 20. The channel filter 22 is a band pass filter whose pass range 

is a frequency band of ±10 MHz around the intermediate frequency. 

The output of the band pass filter is connected to a second mixer 
stage having a second and a third mixer 24 and 26 respectively, to which 
the output signal of the band pass filter is passed in parallel. The two 

20 mixers 24 and 26 are connected at their second input to a second local 

oscillator 28 which is tuned to the intermediate frequency fi F - Provided 
between the output of the local oscillator 28 and the third mixer 26 is a 
phase shifter 30 which causes a phase shift of the oscillation of the local 
oscillator 20 through 90°. 

25 The second mixer stage on the one hand provides for down- 

conversion of the signal, whereby the signal is freed of the frequency 
components of the intermediate frequency. On the other hand, the 
phase-shifted mixing in the parallel mixers 24 and 26 provides for 
separation of the real and imaginary parts of the signal, whereby it is 

30 possible to ascertain the phase position of the signal oscillation. 

Low pass filters (LPA) 32 and 34 connected on the output sides of 
the mixers 24 and 26 provide for anti-aliasing, that is to say smoothing of 



-5- 



the signal components. Finally there are provided analog-to-digital 
converters (ADC) 36 and 38 which are clocked in parallel by a sampling 
oscillator or sampling clock 40. The sampling frequency is f s . The 
accuracy of the sampling frequency is ±f S£ . 
5 Figure 2 shows a diagram illustrating the situation of erroneous 

estimation of the carrier frequency and phase noise. The frequency f of 
the oscillation of the first local oscillator 20 is plotted on an abscissa 50. 
The amplitude of the oscillation is plotted on an ordinate 52. The ideal 
frequency spectrum of the local oscillator (LO) is represented at an 

10 abscissa position f L0 . It is in the form of a delta function, as is illustrated 
here in the form of an arrow 54. In other words, the ideal frequency of 
the local oscillator is exactly the difference of the carrier frequency f c and 
the intermediate frequency fi F . 

A line shape 56 whose maximum is to found on the abscissa at the 

15 position f LO + f e represents the real frequency spectrum of the local 
oscillator 20. The maximum of the line shape 56 is shifted relative to the 
ideal frequency spectrum, by erroneous estimation of the carrier 
frequency. The line shape is further broadened by phase noise. 

Figure 3 shows by way of three symbols selected by way of 

20 example the residual phase error O error which, subsequently to 
synchronization and Fourier transform of a symbol outputted by the down- 
converter in Figure 1, remains in its partial signals. The diagram in Figure 
3 represents the dependency of the residual phase error on the ordinal 
number of the partial signals, by reference to three examples. To explain 

25 the abscissa division: each symbol of a frame is transmitted in the OFDM 
system in the form of partial signals on the predetermined number of 
subcarriers. Associated with each subcarrier i is a predetermined 
frequency interval with a center frequency fj which is clearly determined at 
a given spacing 1/AT of the center frequencies by specifying an ordinal 

30 number i. The abscissa division in Figure 2 is not shown exactly true to 
scale. In the diagram in Figure 3 ordinal numbers from -26 to +26 
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correspond to a total number of 52 partial signals. The ordinal number 0 
is not allocated. 

The phase error can be determined on the basis of 
predetermination by the relevant standards only by means of 
5 predetermined pilot partial signals, referred hereinafter as pilot signals. 

Previously known information which is not related to the useful 
information is transmitted as a pilot signal. The example in Figure 3 has 
four mutually spaced pilot signals distributed over a symbol. The pilot 
signals have in one symbol the ordinal numbers -21, -7, +7 and +21. 
10 The underlying analysis of the present invention affords the 

following main error sources for phase errors in the Fourier-transformed 
signal: 

1. Phase variations due to an erroneous estimation of the phase 
timings. 

15 2. Phase variations due to an erroneous estimation of the 

carrier frequency. 

3. Phase variations due to phase noise. 

4. Phase variations due to a frequency offset of the sampling 
oscillator frequency during analog-to-digital conversion. 

20 The specified sources of phase errors are described in greater detail 

hereinafter. 

1. Erroneous estimation of the frame timing 

In the synchronization operation the first sample signal or sample 

of a frame must be estimated in regard to correct decoding of the data. 
25 In the case of sampled digital data the resolution when determining the 

sample is greatly dependent on the sampling interval. The sampling 

interval is the inverse of the sampling frequency. The timing error can be 

a multiple of the timing interval. The uncertainty in terms of estimation is 

therefore ± 0.5 T S/ wherein T s is the sampling time. 
30 A timing error will appear on the receiver side after performing the 

FFT as a linear phase error. That timing error is the same for all OFDM 
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symbols in the same frame. Therefore all symbols of a frame are 
influenced by a linear phase component in the same manner: 

<D te (k) = ±nk/N, (1) 
In that equation k denotes the frequency at the output of the FFT in the 
5 digital domain, <X> tE (k) denotes the phase error and N denotes the total 

number of the subcarriers. 

2. Erroneous estimation of the carrier frequency 

Upon reverse conversion of the received radio frequency signals the 
local oscillator (LO) is not exactly tuned to the expected frequency. That 

10 frequency offset can be estimated by means of preamble symbols 
provided in each frame. As however that operation is performed in the 
digital domain the estimation is not accurate. Besides thermal noise, 
digital noise which adversely affects the estimation also occurs because of 
the limited number of bits which represent the preambled symbols and 

15 which are used for representation of the frequency offset. 

The frequency offset is usually specified in standardized form with 
respect to the frequency spacing in the OFDM system. In the case of the 
above-mentioned standards the channel spacing is for example Af = 312.5 
KHz. A frequency offset f E of 200 KHz is accordingly specified as a 

20 frequency offset of 0.64. 

A good synchronization unit estimates the frequency offset with an 
error of ±0.001 (standardized). That slight residual frequency error 
occurs in the OFDM symbol as a constant phase after implementation of 
the FFT. The error accumulates from one symbol to another and therefore 

25 produces a large phase shift after a number of symbols. 

3. Phase variations due to phase noise 

The signal generated by the local oscillator is broadened in the 
frequency space due to the phase noise. 

Essentially two effects are linked to the phase noise. The first is the 
30 interference between the carriers (inter-carrier interference or ICI). The 

second effect is what is referred to as the common phase error or CPE. 
ICI appears as additionally Gaussian noise. This can only be reduced by 
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the local oscillator being designed differently. CPE produces a phase shift 
in the OFDM symbol after the FFT. That effect is similar to that of a 
frequency offset. The constant phase shift however is not accumulated 
here. It varies randomly from one symbol to another. 
5 ICI and CPE both depend on the number of subcarriers which are 

used in the OFDM system. In the case of the ICI the power increases 
with the number of subcarriers. The ICI is therefore the limiting factor in 
systems such as DVB or DAB where thousands of subcarriers are used. In 
contrast thereto the CPE increases with a decreasing number of 

10 subcarriers. Therefore CPE represents the limiting factor in systems in 
accordance with the standard IEEE 802.11a and Hiperlan/2 where only 64 
subcarriers are used. The correlation of the phase shifts over various 
symbols is very slight in the case of the CPE. Therefore there is no linear 
prediction method. 

15 4. Phase shifts due to an error in the sampling frequency 

At the receiver input signals are sampled at a given frequency and 
transferred into digital signals. In a real system the actual sampling 
frequency does not exactly conform with the expected frequency. Rather, 
it will have an error f Se which is measured in parts per million (ppm) of the 

20 original sampling frequency. With a sampling rate of 20 MHz with a 

frequency error of ±20 ppm the actual sampling frequency will therefore 
be in the range of 20 MHz ± 400 Hz. 

The four error sources in respect of the phase position referred to 
above afford a linear dependency of the phase <fo(f) on the frequency f in 

25 the form 

4>,(f) = mff + c, (2) 
In that equation the letter i denotes the ordinal number of a symbol, mj 
and q are constants, mj being the gradient and q being a frequency- 
independent phase offset. 
30 The gradient mi is dependent on the ordinal number of the 

respective symbol and can be represented as: 

m, = rn 0 + i£ (3) 
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In that equation £ denotes the component which is caused by the error in 
respect of the sampling frequency. £ can be positive or negative, 
depending on whether the actual sampling frequency is greater than or 
less than the expected value. That error component increases with the 
ordinal number i of the respective symbol. It therefore accumulates from 
one symbol to another. 

m 0 denotes the component which is produced by the error in the 
frame timing. It is the same for all symbols. 

The coefficient q can be represented as 

q = ic 0 + oij (4) 

In that equation otj is the component of the common phase error CPE for 
the symbol in question and c 0 is the phase error which is caused by the 
frequency offset of the carrier frequency. That error accumulates from 
one symbol to another. 

Analysis shows that the sampling frequency error, due to the 
accumulation from one symbol to another, obtains a particularly great 
weight in the overall phase error. As it changes with each symbol, 
correction thereof is particularly difficult. In the method according to the 
invention therefore the equalization step includes a step of at least 
partially eliminating an accumulation of a phase error caused by a 
sampling frequency error in the partial signal over the sequence of the 
partial signals in such a way that the accumulation is negligible. The 
phase error caused by the sampling frequency can be disregarded by 
virtue of that measure in terms of correction of the remaining phase error. 

In the method according to the invention the estimation step 
includes a step of detecting a plurality of predetermined pilot signals and a 
step of determining a phase correction factor on the basis of the detected 
pilot signals, wherein at least one multiplication operation is performed 
just by means of shift and adding operations. 

Shift and adding operations are particularly simple operations with 
digital data, which can be performed at great speed. The idea of the 
invention of using shift and adding operations in determining the phase 
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correction factor is closely linked in the method according to the invention 
to the previously stated method steps. The possibility of disregarding the 
phase error caused by the sampling frequency affords the prerequisite for 
simplifying estimation of the phase error in such a way that its 
5 mathematical calculation can be reduced at least predominantly and, in a 
particularly preferred embodiment described hereinafter, even completely, 
to the implementation of shift and adding operations, using the previously 
determined pilot signals. 

Accordingly the method according to the invention achieves a high 

10 speed in determining the phase error, in comparison with known methods. 

Just the replacement of a multiplication operation by shift and adding 
operations affords speed advantages over known methods. In that way, 
the delay in the spread of the received signal, which is caused by the use 
of such a method, becomes particularly slight. That is of significance in 

15 particular in communication processes by way of wireless channels. 

It is a basic point that the speed of the method according to the 
invention increases, the fewer real multiplication operations are used. 
Therefore, in accordance with the invention, the optimum in terms of 
speed is achieved if no multiplication operations are involved. When lower 

20 requirements in respect of speed are made however, it is also possible to 
use multiplication operations within the context of the invention. 

The method according to the invention is particularly suitable for 
use in accordance with the standards IEEE 802.11 and Hiperlan/2. In that 
respect it is assumed in particular that the subcarriers are not subjected 

25 to great changes from one symbol to another, as is the case with a 
maximum speed of mobile transmitters and receivers of 3 m/s. 

Various embodiments of the method according to the invention are 
possible, in which the phase error, due to the sampling frequency 
deviation, can be rendered negligible. 

30 In a preferred embodiment of the method according to the 

invention the equalization step includes a step of division of the partial 
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signal by a complex second signal which has the phase of the partial 
signal preceding in the sequence of the partial signals. 

The above analysis affords the following for the phase shift of a 
symbol i: 

<)>i(k) = (m 0 + i£)k + ic 0 + (Xj (5) 
and the symbol i+1 

+i + i(k) = (m 0 + (i+l)^)k + (i+l)c 0 + <x, +1 (6) 
It will be seen that phase error r£ caused by the sampling error 
rises with an increasing ordinal number i in the sequence of the partial 
signals. The phases are subtracted by complex division. That therefore 
gives the following: 

+i+i(k) - <h(k) = S'k + Co + (<x, +1 -a.) (7) 
The difference in the phase shifts is therefore not dependent on the 
ordinal number i. That means that the accumulation of the phase error, 
due to the sampling frequency, can be eliminated by means of division of 
the current complex signal by a second complex signal if that second 
complex signal contains the phase information of the preceding signal in 
the sequence. 

The operation of ascertaining the second signal is the content of 
channel estimation. Preferably the second signal is the frequency 
response function of the associated subcarrier. The frequency response 
function of the subcarrier contains information about the attenuation and 
the phase shift of the subcarrier. In that way, in the equalization step, it 
is possible to effect a phase correction which makes it possible to 
substantially eliminate the phase error caused by erroneous estimation of 
the sampling frequency. 

In this embodiment, there is preferably provided for each partial 
signal with the exception of the pilot partial signals a step of determining 
the frequency response function on the basis of a preceding partial signal 
in the sequence of the partial signals. The cyclically repeated operation of 
determining the frequency response function is based on the realization 
that the accumulation of the phase error due to the sampling frequency 
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deviation, can be best prevented by determining the current phase shift 
on the channel. Basically it is also conceivable for the frequency response 
to be determined not for each partial signal but for example only for each 
second partial signal, and for an interpolated signal to be used in the 
5 meantime. In that way however the accumulation is not completely 

eliminated and the quality of phase correction is worsened. The 
consequence of this is that the error rate in decoding is increased. 

The step of determining the frequency response function preferably 
includes a division of the preceding partial signal in the sequence by a 
10 third signal, wherein the third signal represents the equalized, phase- 

corrected and decoded and then recoded preceding partial signal in the 
sequence. The third signal thus represents the preceding partial signal 
cleared of all errors. After the division operation attenuation of the 
subcarrier can be inferred from the value of the quotient and the phase 
15 shift of the subcarrier can be inferred from the phase of the quotient. The 
phase shift can be determined by a procedure whereby the arc tan of the 
quotient of the imaginary part and the real part is subtracted from the 
phase of the quotient. 

Preferably the method according to the invention provides 52 
20 subcarriers. That corresponds to the number of subcarriers in the 
standards IEEE 802.11a) and Hiperlan/2. 

In this respect there are preferably provided four pilot subcarriers 
also corresponding to the stated standards. If the numbering of the 
subcarriers: -26, -25, -24, -1, 1, 2, 3, 26 is followed, then the 
25 ordinal numbers -21, -7, 7 and 21 are associated with the pilot subcarriers 

in this embodiment. 

In a particularly preferred embodiment of the invention the 
estimation step includes a step of calculating the following parameters on 
the basis of the pilot signals: 
30 cos(<M = ( 1 /4)(tte{P- 2l } + *e{P_ 7 } + *e{P +7 } + *e{P +21 }), (8) 

sin((|)o) = ( 1 A)(3m{P. 2l } + 3m{P. 7 } + 3m{P +7 } + 3m{P +21 }) / (9) 
27ip 0 sin((f>o) m (2^e{P.2i>+3^e{P. 7 }-3lRe{P + 7}-2lHe{P + 2i})/128 / (10) 
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-27rp 0 cos((t)o) « (23m{P. 2 i}+33m{P- 7 }-33m{P +7 }-23m{P +21 })/128 / (11) 
Therein P. 2 i, P-?, P+?/ P+21 denote pilot signals and <Ke and 3m denote the 
operation of determining the real and imaginary part respectively. 

Calculation of those parameters is based on the following model: 
The signal of the pilot channels can be represented as follows: 

p k = e j(2*pok + * 0 ) ; j = ^/n (12) 

Therein p 0 denotes the frequency offset of the sampling oscillator 
(sampling clock frequency offset) and is directly proportional to the 
above-discussed factor $ 0 contains the common phase error due to 
phase noise and any residual error in the carrier frequency estimation. In 
equation (12) amplitude changes in the channel are disregarded. Only 
phase changes are taken into consideration. 

On the basis of equation (12) once again to compare the solution 
according to the invention to the state of the art it is pointed out that, in 
the state of the art, the direct calculation of <J> 0 and p 0 is effected by 
resolving equation (12) in accordance with those parameters. That 
requires a high level of computing complication and expenditure, in 
particular arc tan calculation and correction by means of an NCO block. 
Both operations can be implemented in a CORDIC algorithm. 
Under the condition: 

27ikp 0 <<1 (13) 
P k can in a first approximation be approximated to the exponential 
function as follows: 

P k = cos(<t> 0 ) - 27ip 0 k sin(<|>o) + j(sin(4>o) + 27ip 0 k cos((|>o)) (14) 

That affords the following for the pilot signals: 
P-21 = cos(<|>o) + 27ip 0 21 sin(4)o) + j(sin(<J) 0 ) - 27ip 0 21 cos(<t> 0 )) (15) 
P-7 = cos(<|>o) + 27cp 0 7 sin(<J)o) + j(sin(<t> 0 ) - 2np 0 7 cos{M) (16) 
P+7 = cos(<|>o) - 27ip 0 7 sin(<|)o) + j(sin(<|)o) + 2np 0 7 cos(^ 0 )) (17) 
P+21 = cos(<|)o) - 27ip 0 21 sin(<j> 0 ) + j(sin(<t> 0 ) + 27ip 0 21 cos((|)o)) (18) 
Accordingly the following four parameters have to be determined: 
cos(4> 0 ), sin^o), 2rcp 0 sin((|>o), 27ip 0 cos((|)o). That is effected for the first two 
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10 



15 



20 



25 



parameters without multiplication solely on the basis of shift and adding 
operations in the following fashion: 

cos(<|>o) = (V4)(5Re-CP. 21 } + SRe{P. 7 } + tte{P +7 } + tte{P +21 }) (19) 
sin(4> 0 ) = ( 1 /4)(3m{P. 2l } + Dm{P. 7 } + 3m{P +7 } + 3m{P +21 }) (20) 
A further approximation is introduced for the last two parameters, 
to avoid multiplication operations. Combining the real parts of the pilot 
signals firstly gives the following: 
2tc Po 35 sin(4> 0 ) = (V 2 )(9le<P. 21 }+2$Re{P. 7 }-29?e{P +7 }-5Re{P +21 }) (21) 
27tp 0 28 sin((M = ( 1 /2)(9le{P. 2l }+9le-CP- 7 }-lRe-CP +7 }-9te{P +2l }) (22) 
In equation (21) the number 35 is now approximated by the 
number 36. Using the following two reductions: 

4^ 



36 = ( 32 + 
32 t 32 



and 



28 
32 



32-4 
32 



V 



8 



(23) 



(24) 



that gives (25) 

2kp 0 sin(<>o)|| + 2*p 0 sin(<t> 0 )|| * 2np 0 sin(* 0 )|| + 2tup 0 sin(<|>o)|| = 2np 0 sinfoo) • ^ 



1 + 



i) 



and finally 

27cp 0 sin((|)o) « (2«e{P.2i}+3lRe{P.7}-39?e{P +7 }-2^e{P +21 })/128 (26) 
The exact value of the foregoing expression is 

(126/128)'27ip 0 sin((t>o). Accordingly the approximation involved gives rise 

to an error of about 1.6%. 

The following is obtained in a similar manner: 

-27ip 0 cos(4>o) « (23m{P. 21 }+33m{P. 7 }-33m{P +7 }-23m{P +21 })/128 (27) 

Therefore, with an only slight error, the present embodiment makes 

it possible to totally avoid multiplication operations for determining the 

parameters. 

In a further preferred method the estimation step includes a step 
of storing the parameters 2np 0 s\r\($ 0 ) and (-27tp 0 cos((|)o)) in a respective 
first register. 
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In this embodiment of the invention the estimation step includes a 
first step of multiplying the respective content of the first register by a 
factor 26, that multiplication operation being carried out by means of a 
reduction 26 = 2 5 - 2 2 - 2 in the form of shift and adding operations. 
5 That method is based on the following correction algorithm: 

After detection of the pilot signals, partial signals with useful data 
are detected. The first partial signal is identified by D_ 26 . 

The phase component on the subcarrier in question, due to the 
residual phase error, can be represented as follows: 
10 D-26 = cos(c|)o) + 27ip 0 26 sin(<|)o) + j(sin(<|>o) - 27ip 0 26 cos(()>o)) (28) 

The required correction is accordingly D*- 26 = conj{D. 26 }, wherein 
conj{.} denotes the complexly conjugated partial signal. Two 
multiplications would be required in themselves in order to determine that 
phase correction factor in accordance with equation 28. By using the 
15 reduction 26 = 2 5 - 2 2 - 2 however the multiplication operations can be 

replaced by shift and adding operations. Preferably the product calculated 
in that way is written into a second register. 

In this embodiment, in a particularly preferred configuration, the 
phase correction factors of the following partial signals can also be 
20 determined without multiplication operations. In the case of a following 
partial signal the correction step, after the above-described step of 
"multiplication", includes a step of inverting the content of the first 
register, a step of adding the inverted content of the first register and the 
content of the second register, and a step of overwriting the second 
25 register with the calculated sum. 

This is explained by the example of the correction factor D*_ 25 . The 
following occur here: 

<Ke{D*. 25 } = lHe{D*. 26 }-27ipoSin((|)o) (29. a) 

3m{D*- 25 } = 3m{D*. 26 }-(-27cp 0 cos((t)o)) (29. b) 

30 Accordingly D*. 24 arises out of D*_ 25 and so forth. It is apparent that 

in this way the phase correction factors can be determined from the 
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respectively preceding ones without an actual multiplication operation 
having to be performed. 

Particular attention is required in this embodiment in respect of the 
phase correction factors which directly follow a pilot signal from the 
5 ordinal number: D*. 2 o, D*. 6 , D* +1 , D* +8 and D*+ 2 2. In those cases a 

"sampling hole" is given by a partial signal. Those phase correction 
factors are therefore determined as follows: 





*e{D*. 20 } = 


: 9te{D*- 22 }-27tp 0 2sin(<t> 0 ) 


(30.a) 




3m{D*. 20 } = 


= 3m{D*- 22 }-(-27tp 0 2cos(<|>o)) 


(30. b) 


10 


<Re{D*_ 6 } = 


9?e-CD*. 8 }-27tp 0 2sin(<t>o) 


(30.c) 




3m{D*. 6 } = 


3m{D*- 8 }-(-27tpo2cos(<|>o)) 


(30.d) 




<Re{D* +1 } = 


5Re{D*.i}-27cp 0 2sin((()o) 


(30.e) 




3m{D* +1 } = 


= 3m{D*. 1 }-(-27tp 0 2cos(<|)o)) 


(30.f) 




*e{D\ 8 } = 


9?e{D* +6 }-2jip 0 2sin((t)o) 


(30.g) 


15 


3m{D\ 8 } = 


= 3m{D* +6 }-(-27tp 0 2cos(<t)o)) 


(30. h) 




*e{D* +22 } = 


= 9?e{D* +20 }-27ip 0 2sin(<j)o) 


(30.i) 




3m{D* +22 } 


= 3m{D* +20 }-(-2np 0 2cos((()o)) 


(30.j) 



In regard to the apparatus aspect the technical object of the 
invention is attained by a phase correction unit including a first signal 
20 input and a computing unit for determining a complex phase correction 
factor by means of arithmetic operations which include multiplication. 
According to the invention there is provided a second signal input to which 
the computing unit is connected, and in addition at least one multiplication 
operation is implemented in the form of at least one shift operator in 
25 conjunction with at least one adder. 

In the case of the phase correction unit according to the invention, 
there is provided a second signal input for detecting pilot signals. 

The phase correction unit according to the invention permits fast 
correction of the phase of a partial signal, as was described hereinbefore 
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in connection with the method of the invention. In regard to the 
advantages of the apparatus therefore attention is expressly directed to 
the description of the method aspect of the invention. 

Preferably all multiplication operations are implemented in the form 
5 of at least one shift operator in conjunction with at least one adder. 
Multiplication with an invariable factor is preferably implemented in the 
form of a reduction of the factor in powers of 2 of the kind a 0 2° + ai 2 1 + 
a 2 2 2 + with a 0/ ai, a 2/ ... = 1 or -1. In that respect the example of 
multiplication with a factor 26 has been described hereinbefore. 

10 The computing unit is preferably adapted to determine the 

parameters of equations 21, 22, 26 and 27 with the signals (P- 2 i/ P-7, P+7, 
P+21) received by way of the second signal input. In that case there is 
preferably provided a specific parameter calculation unit of suitable 
configuration, for determining each of the parameters. 

15 The invention is described hereinafter by means of an embodiment 

by way of example. For that purpose reference is made to Figures 4 
through 7 which are firstly briefly summarized jointly with above- 
described Figures 1 through 3. In the drawing: 

Figure 1 shows a simplified block circuit diagram of a down- 

20 converter in accordance with the state of the art, 

Figure 2 shows a diagram to illustrate a frequency offset and the 
effect of phase noise, 

Figure 3 shows a diagram to illustrate the accumulation of the 
phase error in symbols after the Fourier transform, 

25 Figure 4 shows a block diagram of an embodiment of an OFDM 

receiver, 

Figure 5 shows a detailed block diagram of the OFDM receiver of 
Figure 4, 

Figure 6 shows a detailed block diagram of the phase correction 
30 unit of Figure 5, and 

Figure 7 shows a detailed block diagram of the accumulator of 
Figure 6. 
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Figure 4 shows a block diagram of an embodiment by way of 
example of an OFDM receiver 70. The OFDM receiver 70 has a unit 72 for 
extracting cyclic prefix signals from the received signal. Connected on the 
output side thereof is an FFT block 74 which performs a fast Fourier 
5 transform. The output of the FFT block 74 is connected on the one hand 
to an equalization block 76 which is also referred to hereinafter as the 
equalizer. On the other hand the output of the FFT block 74 is connected 
to a channel estimation block 78 whose output is in turn taken to the 
equalizer 76. 

10 The output of the equalizer is connected to a phase correction block 

80 which in turn communicates with a block 82 for demapping and for 
forward error correction (FEC). 

Figure 5 shows a detailed block diagram of the OFDM receiver 70. 
Fewer details than in Figure 4 are shown only at the input of the OFDM 

15 receiver 70. Here the unit 72 is combined together with a synchronizer 

and the FFT block in a single block 84. 

The block 82 in Figure 4 includes per se known components: a soft 
demapper 82, a deinterleaver 88 and a soft Viterbi decoder 90. The 
output of the soft Viterbi decoder 90 is fed back to an encoder 92, on the 

20 output side of which are connected an interleaver 94 and a mapper 96. 

A zero-forcing estimator 98 is connected to the mapper 96. 
Connected on the output side thereof is a buffer memory 100. The buffer 
memory 100 has a writing control 102, an S-RAM memory 104 and a 
reading control 106. 

25 Provided between the output of the FFT block 74 and the writing 

control 102 is a first sign unit 108 for correcting the sign of reference 
signals. Also provided between the output of the FFT block and the zero- 
forcing estimator 98 is a buffer memory in order to produce a defined 
delay of data which are sent from the FFT block 74 to the zero-forcing 

30 estimator 98. Finally provided between the buffer memory 110 and the 
writing control 102 is a second sign unit 112 which sets the sign of pilot 
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signals positive, independent of the sign with which they occur at the 
input. 

The mode of operation of the OFDM receiver 70 is described in 
greater detail hereinafter. Firstly channel estimation in the block 78 is 
5 described in greater detail. In this case, the standards IEEE 802.11a and 
Hiperlan/2 form the basis in terms of the structure of the received signals. 

The first symbol received from the synchronizer 84 is the reference 
symbol H REF (k). The reference symbol is derived from the long preamble 
symbols. That symbol is used by the receiver as a reference because the 
10 data communicated therein are known to the receiver from outset. 

The reference symbol H REF (k) is then stored in the buffer memory 
100. In that respect k denotes the subcarrier frequency with the following 
values (in this order:) k = -21, -7, +7, +21, -26, -25, -22, -20, -19, 
-9, -8, -6, -5, -1, +1, +5, +6, +8, +9, +19, +20, +22, 
15 +25, +26. The order of the subcarriers is produced in the FFT block 74. 

In all accordingly 52 signals are stored. 

The first symbol after the reference symbol, the so-called SIGNAL 
symbol, is then corrected by the equalizer 76. The equalizer 76 calculates 
the quotient of the SIGNAL symbol and the symbol stored in the buffer 
20 memory 100. 

The symbol equalized by the equalizer 76 is then subjected to 
phase correction in the phase correction block 80, completely 
demodulated and finally decoded. 

The bits which are decoded from the SIGNAL symbol in the Viterbi 
25 decoder 90 provide information about the rate and the number of bytes 
which are received in the current frame. The SIGNAL symbol is not read 
into the feedback loop beginning with the encoder 92. 

The first DATA symbol (symbol #1) which is received after the 
SIGNAL symbol is then corrected again in the equalizer 76 with the 
30 channel estimation contained in the buffer memory 100. 

The symbol outputted by the equalizer is phase-corrected, 
demodulated and then decoded. The bits at the output of the Viterbi 
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decoder, which correspond to the symbol #1 are fed back into the channel 
estimation block 78. 

In the channel estimation block 78 the bits are firstly recoded and 
then modulated. The complex signals X(i-l,k) obtained are compared to 
5 the complex signals which were present for the same symbol at the stage 
prior to equalization in the equalizer. This comparison involves a division 
operation. That form of channel estimation therefore represents a zero- 
forcing estimation. That channel estimation is stored in the buffer 
memory 100. 

10 The above-specified division is performed for 48 partial signals, 

namely those which correspond to k = -26, -25, -22, -20, -19, -9, - 
8, -6, -5, -1, +1, +5, +6, +8, +9, +19, +20, +22, +25, +26. 
The partial signals for k = -21, -7, +7 and +21 (pilot signals) are stored 
directly in the buffer memory 100 after sign correction in the block 112. 

15 The channel estimation obtained in that way acquires the 

information about the residual phase error because the data in the 
feedback branch do not have any phase error at all. 

The next received symbol is treated precisely like the preceding 
one, beginning with equalization in the equalizer 76. 

20 Reference is now made to Figure 6 to describe in greater detail 

phase correction in the phase correction block 80. For that purpose the 
structure of the phase correction block 80 will firstly be described. 

At the input side the phase correction block 80 has a switch 120 
controlled by a switch control 122. The position of the switch, which is 

25 the "upper" position in Figure 6, provides for a conducting connection to a 
unit 124 which forwards the real and the imaginary parts of the inputted 
signal separately into two parallel processing branches which are 
described hereinafter. 

The real part is passed to two parameter circuits 126 and 130. The 

30 first parameter circuit 126 determines the parameter cos O 0 in accordance 
with equation 8. The second parameter circuit 130 determines the 
parameter in accordance with 27ip 0 sin O 0 in accordance with equation 10. 
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The imaginary part is passed in parallel to a third and a fourth parameter 
circuit 128 and 132 respectively. The third parameter circuit determines 
the parameter sin O 0 in accordance with equation 9. The fourth 
parameter circuit 132 determines the parameter -27cp 0 sin O 0 . 
5 A respective accumulator 134 and 136 is connected on the output 

side of the second and fourth parameter circuits 130 and 132. The 
structure and function of the two accumulators which are of an identical 
design are described hereinafter with reference to Figure 7. 

The outputs of the third parameter circuit 128 and of the 

10 accumulator 136 are fed to an adder 138. The outputs of the first 
parameter circuit 126 and the accumulator 134 are fed to an adder 140. 

The output of the adder 138 is connected by way of an inverter 142 
to an imaginary part input of a circuit 144 whose real part input is 
connected to the output of the adder 140. The inverter 142 inverts the 

15 sign of the output signal of the adder 138. The circuit 144 converts the 
signals at the two inputs into a complex signal which, like the signal at the 
second output of the switch 120, is passed to a complex multiplier 146. 

The structure of the accumulators 134 and 136 is described 
hereinafter with reference to Figure 7. Figure 7 admittedly relates to the 

20 accumulator 134 but the accumulator 136 is of the same structure. 

The accumulator 134 has a first register 150 into which signals are 
written from the parameter circuit 130. The output of the register 150 is 
passed in parallel on the one hand to a factor circuit 156 and on the other 
hand to a second register 154. The factor circuit 156 multiplies the output 

25 signal of the first register by 26. That is effected solely by means of shift 
and adding operations on the basis of the above-indicated reduction of the 
number 26 into powers of the number 2. The output of the factor circuit 
156 is connected to a third register 152 by way of a first input of a switch 
158 controlled by a control circuit 160. 

30 The output of the second register 154 is passed in parallel to two 

inputs of a switch 162, in which case however one of the two parallel 
branches is passed by way of a factor circuit 164 which by means of a 
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shift operation effects multiplication of the register output signal by 2. 
The output of the switch 162 and the output of the register 152 are 
connected to the inputs of an adder 166. 

The circuit in Figure 6 calculates for a respective signal the phase 
correction factor, as for example corresponding to equations 29. a and 
29. b. The subtraction operation required for the respective symbol is 
effected by means of the accumulator 166. As discussed above the value 
to be subtracted is dependent on the ordinal number of the symbol to be 
corrected. 

The circuit, by suitable control of the switches in the accumulator 
circuit in Figure 7, also permits calculation of the "special cases" 
corresponding to equations 30. a through 30.j. 

That circuit affords a particularly simple structure for a phase 
correction circuit which causes a particularly slight delay. 
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CLAIMS 



1. (currently amended) A method of reducing a phase error 
caused by a plurality of error sources in a signal which is present in a 
digital frequency representation in the form of a sequence of a plurality of 
digital partial signals which are associated with a number of subcarriers 
(k) of a carrier, with the steps which or e effected comprising, for each 
partial signal , the steps of : 

- equalization of the partial signal (Y(i,k)), 

- estimation of the phase error of the equalized partial signal 
(X(i,k)), and 

- correction of the estimated phase error of the equalized partial 

signal, 

characterized in thati 

- the equalization step includes a step of at least partially 
eliminating an accumulation of a phase error of the partial signal, caused 
by a sampling frequency error, over the sequence of the partial signals, 
such that the accumulation is negligible, and 

- the estimation step includes a step of detecting a plurality of 
predetermined pilot signals and a step of determining a phase correction 
factor on the basis of the detected pilot signals, wherein at least one 
multiplication operation is carried out solely by means of shift and adding 
operations. 

2. (original) A method as set forth in claim 1 characterized in 
that the equalization step includes a step of division of the partial signal 
by a complex second signal which has the phase of the partial signal 
preceding in the sequence of the partial signals. 

3. (original) A method as set forth in claim 2 characterized in 
that the second signal represents a frequency response function of the 
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associated subcarrier for the preceding partial signal in the sequence of 
partial signals. 

4. (original) A method as set forth in claim 3 characterized by a 
step, effected for each partial signal with the exception of the pilot partial 
signals, of determining the frequency response function on the basis of a 
preceding partial signal in the sequence of partial signals. 

5. (original) A method as set forth in claim 4 characterized in 
that the step of determining the frequency response function includes a 
division of the preceding partial signal in the sequence by a third signal, 
wherein the third signal represents the equalized, phase-corrected and 
decoded and then recoded preceding partial signal. 

6. (currently amended) A method as set forth in one of the 
preceding c l aims claim 1 characterized in that there are provided 52 
subcarriers. 

7. (original) A method as set forth in claim 6 characterized in 
that there are provided four pilot subcarriers for the transmission of pilot 
signals with which the ordinal numbers -21, -7, 7 and 21 are associated 
with the following numbering of the subcarriers: -26, -25, -24, -1, 1, 2, 
3, 26. 

8. (original) A method as set forth in claim 7 characterized in 
that the pilot signals are present at the beginning of the sequence of 
partial signals. 

9. (original) A method as set forth in claim 8 characterized in 
that the phase correction factor is determined utilizing the pilot partial 
signals without a multiplication operation. 

10. (currently amended) A method as set forth in one of the 
preceding c l aims claim 1 characterized in that the correction step includes 
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a step of complex multiplication of the equalized partial signal with the 
phase correction factor. 

11. (currently amended) A method as set forth in c l aim 0 or 
claim 10 characterized in that the estimation step includes a step of 
calculating the following parameters on the basis of the pilot signals: 

cos(<t) 0 ) = ( 1 / 4 )0Re{P- 2l } + *e{P. 7 } + tte{P +7 } + tte{P +21 }), 
sin((t> 0 ) = ( 1 /4)(3m{P. 2l } + 3m{P. 7 } + 3m{P +7 } + 3m{P +21 }), 
27ip 0 sin((|)o) « (29?e{P. 21 } + 39ie{P- 7 }-3lRe{P +7 }-2^e{P +2x })/128 / 
-27ip 0 cos((|>o) * (23m{P. 21 } + 33m{P. 7 }-33m{P +7 }-23m{P +21 })/128, 

wherein P. 2 i, P. 7/ P +7 , P+ 21 denote pilot signals and 9Te and 3m denote the 
operation of determining the real and imaginary part respectively. 

12. (original) A method as set forth in claim 9 characterized in 
that the estimation step includes a step of storing the parameters 
27rp 0 sin((J>o) and -27tp 0 cos((|>o) in a respective first register. 

13. (original) A method as set forth in claim 12 characterized in 
that the estimation step includes a first step of multiplying the respective 
content of the first register by a factor 26, wherein said multiplication 
operation is effected by means of a reduction 26 = 2 5 - 2 2 - 2 in the form 
of shift and adding operations. 

14. (original) A method as set forth in claim 13 characterized in 
that the calculated product is written into a second register. 

15. (original) A method as set forth in claim 14 characterized in 
that in the case of a subsequent partial signal the correction step after the 
multiplication step includes a step of inverting the content of the first 
register, a step of adding the inverted content of the first register and the 
content of the second register, and a step of overwriting the second 
register with the calculated sum. 
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16. (original) A phase correction unit including a first signal input 
and a computing unit for determining a complex phase correction factor 
by means of arithmetic operations which include multiplication, 
characterized in that there is provided a second signal input to which the 
computing unit is connected, and at least one multiplication operation is 
implemented in the form of at least one shift operator in conjunction with 
at least one adder. 

17. (original) A phase correction unit as set forth in claim 16 
characterized in that all multiplication operations are implemented in the 
form of at least one shift operator in conjunction with at least one adder. 

18. (currently amended) A phase correction unit as set forth in 
claim 16 or cloim 17 characterized in that multiplication with an invariable 
factor is implemented in the form of a reduction of the factor in powers of 
2 of the kind a 0 2° + a x 2 1 + a 2 2 2 + with a 0/ ai, a 2/ ... = 1 or -1. 

19. (currently amended) A phase correction unit as set forth in 
claim 16 , claim 17 or claim 18 characterized in that the computing unit is 
adapted to determine the following parameters with the signals (P. 2 i, P-7, 
P+7, P+21) received by way of the second signal input: 

cos(4> 0 ) = ( 1 /4)(tte{P- 2l } + *e{P_ 7 } + <He{P+ 7 } + *e{P+ 21 }), 
sin(4>o) = (V4)(3m{P, 21 } + 3m{P_ 7 } + 3m{P +7 } 4- 3m{P+ 21 }), 
27ip 0 sin((t)o) « (29?e{P- 21 }+3^e{P. 7 }-39ie{P + 7}-2^e{P +21 »/128 / 
-27ip 0 cos(<t)o) » (23m{P- 21 }+33m{P. 7 }-33m{P + 7}-23m{P +21 })/128 / 

wherein P- 2i , P. 7/ P+7, P+21 denote pilot signals and <Ke and 3m denote the 
operation of determining the real and imaginary part respectively. 

20. (original) A phase correction unit as set forth in claim 19 
characterized in that for determining each of the parameters there is 
provided a specific parameter calculation circuit of corresponding 
configuration. 
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Abstract 

A method of reducing a phase error caused by a plurality of error sources 
in a signal wh i ch is present i n a digita l frequency representation in the 
form of a sequence of a plurality of digital partial signals which arc 
associated with a number of subcarriers (k) of a carrier- — The fo l lowing 
st e ps arc performed , the method including, for each partial signal: 
equalization of the partial signal (Y(i,k)), estimation of the phase error of 
the equalized partial signal (X(i,k)), and correction of the estimated phase 
error of the equalized partial signal, An- One embodiment of that method 
provides that the equalization step inc l udes with the elimination of an 
accumulation of a phase error of the port i a l signal, couscd by o samp l ing 
frequency error, over the sequence of the partial signals , such that the 
accumulation is neg l igib l e . In addition the estimation step includes a step 
of— detecting a plurality of predetermined pilot signals and a step of 
determining a phase correction factor on the basis of the detected pilot 
signals, wherein with at least one multiplication operation is-carried out 
solely by means of shift and adding operations. A corresponding 
apparatus is also described. 
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